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Abstract 

The simplest anharmonic characteristics of Ir and Rh are discussed in the framework of 
a previously developed simple pseudopotential model which describes the elastic moduli, 
phonon spectra and the lattice heat capacity in the harmonic approximation of these 
metals succesfully. The microscopic Gruneisen parameters, the dependences of the elastic 
moduli on pressure, the coefficient of thermal expansion and the equations of state at the 
finite temperatures have been calculated. The ab initio calculations of the energy-band 
structure and the equation of state for Ir at T = have been done to test the model for 
adequacy at high pressures. The values of different contributions (zero-point oscillations, 
quasiharmonic, etc.) in the considered thermodynamic characteristics of Ir and Rh are 
discussed. 

§1. Introduction 

As it is stated in the literature (Gornostyrev et al 1994, Ivanov et al 1994, 
Katsnelson et al 1996, hereinafter referred to as I), Ir and its analogue Rh 
are set off from other FCC metals because of some distinctive features of 
their physical properties: unusual deformation-induced failure, peculiar de- 
fect structure, specific temperature dependence of the effective Debye tem- 
perature, etc. Moreover, the situation is not trivial because in the case of Ir, 
unlike the most, if not all of the other FCC metals a rather broad range of 
properties (elastic moduli, phonon spectra and lattice heat capacity) can be 
described in terms of a simple pseudopotential model (see I and references 
therein). Even in sp-metals the situation is less favorable: to describe the 
elastic moduli of Ca allowance must be made for the singular contributions 
to energy which result from the proximity of the Fermi surface and the faces 
of the Brillouin zone (Katsnelson et al 1990), whereas for Al it is necessary 
to take into account the contributions of three body forces. Nevertheless, in 
describing the phonon spectra of Al it is impossible to attain an accuracy 
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identical to that for Ir (Browman et al 1971, Fomichev et al 1990). In Ir the 
contribution of three body interactions proves to be small (Greenberg et al 
1990). As a whole, it turns out that surprising as it may seem, Ir is proba- 
bly exceeded in the accuracy of the psedopotential description of the lattice 
properties only by alkali metals. This gives hope that the theory can be 
used to predict those lattice properties of Ir for which the experimental data 
are scanty or absent, such as thermal expansion or the equation of state at 
different temperatures. To be sure, the pseudopotential model is oversimpli- 
fied. Therefore, where it is possible, the model predictions should be tested 
against either the experiment or the results of the ab initio band calculations. 
Despite the fact that no simple model can completely substitute the ab ini- 
tio calculations, to obtain a sufficiently reliable physical information about 
the complex properties of metals even with the help of a model approach 
appears to be of interest. At the same time, fully ab initio calculations of 
the phonon contributions to the thermodynamical properties of metals seem 
for the present to be impossible (or, in any case, are rather cumbersome). 
Therefore, the use of the pseudopotential model for these purposes, where it 
is possible seems to be justified. This paper presents the results of the calcu- 
lations of the coefficient of thermal expansion at finite temperatures, and the 
dependence of the elastic moduli on pressure in terms of the pseudopotential 
model developed previously for Ir and Rh (see I). To understand to what 
extent the simple model used is reliable at high pressures the electron energy 
spectrum of Ir and its equation of state at T = were calculated ab initio . 
The results of the calculations show that the equation of state calculated in 
the framework of the pseudopotential model seems to be sufficiently reliable. 

§2. Thermal expansion and Gruneisen parameters of Ir and Rh 

In calculating the thermal expansion and the Gruneisen parmeters 7(T) 
at finite temperatures as well as the dependencies of the elastic moduli on 
pressure Bik(p) use was made of the pseudopential model constructed in I 
and giving the optimal description of the phonon spectra. The results of the 
calculations in terms of an earlier model (Ivanov et al 1994) are given for 
comparison. The local pseudopotentials V ps (q) used in the two models are 
closely similar. However a slight change in the form of V ps (q) was made in 
I to describe better the phonon spectra. The screening function was taken 
in both cases in the Geldart and Taylor (1970) approximation (see, in more 
detail, Greenberg et al 1990). The temperature dependence of the volume 
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per atom Q(T) = Q + Afi(T) (fi = Q(T = 0)), the coefficient of thermal 
expansion 

and the Gruneisen parameter 7(T) are determined by the following expres- 
sions (Vaks et al 1978): 

Afi(T) _ Pph (T) 
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where 



Cv(T) E(^) *<i + *> 

2V* = 7^ S hu x N x -i\ (5) 

A 

is the phonon pressure; A = qu, q is the phonon wave vector (sweeping the 
Brillouin zone); v is the branch number; uj v is the phonon frequency; 

7A = " m^n (6) 

are the microscopic Gruneisen parameters; 

* = — ETT- (7) 

exp — - 1 

Bt is the bulk modulus at a constant temperature; CV(T) is the phonon 
heat capacity at a constant volume. The microscopic Gruneisen parameters 
are calculated by the formula 

7,(q) = -^I>y*^4, (8) 
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where D is the dynamic matrix; e^ q are the polarization vectors; a and (3 
are the Cartesian indices. In our calculations allowance was also made for 
the quasiharmonic corrections: 



a — a,\ 
Ail _ /Att\ 

7T ~ v7T/i 



1 - 



/AQ\ dlnB 
\~TT)i dlnQ 



/ _8h.il> 



n Ji \ dlntt 



(9) 



where index "1" means the corresponding expression taken in the lowest order 
of anharmonicities. The procedure of calculating the thermal expansion of 
metals and its related quantities was discussed in more detail in Vaks et al 
1978. 

The results of the calculations are presented in figs. 1-4 and table 1. 
First of all, it should be noted once more that 7„(q) vary essentially over 
the Brillouin zone. This indicates again that the Gruneisen approximation 
7„(q) = const is completely unsuited for actual systems. As the temper- 
ature increases the macroscopic Gruneisen parameter 7(T) reaches quickly 
(at T > 0.16, where 6 is the Debye temperature) a constant high temper- 
ature value. Alkali metals exibit the analogous behavior (Vaks et al 1978). 
However, unlike alkali metals having 7 m 1, in Ir and Rh, as in majority of 
transition metals, 7 ~ 2. Note also that the calculated values of 7(T) are 
rather highly sensitive to the form of the pseudopotential used (compare the 
different curves in figs 3b and 4b). 

§3. Equation of states. Used formulas and approximations 

The equation of state for a metal, that is, the dependence of the pressure 
p on Q and T is determined by the expression 

P (n, t) = po(fi) + Pzp {9) + Pph (n, t) + A Pe (n, t) (10) 

where 

E(Q) is the total energy of the crystal at T = 0; 

p ^ = -m^^f) = ^ huj ^ (12) 
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is the pressure relating to the zero-point energy; p ph is the phonon pressure 

(5); 

7T 2 2 dN(Ep) . , 

Ap < = T ~xt (13) 

is the temperature dependent contribution to the electron pressure; N(E F ) 
is the density of states at the Fermi level E F . We shall neglect the anhar- 
monic contributions having the same order of magnitude as the deviation of 
the Ap e (T) from the "low-temperature" (T <C Ep) expression (13). Gener- 
ally speaking, the allowance for the contribution of the zero-point vibrations 
p zp (fl) to the equation of state in the scheme with a pseudopotential whose 
parameters are fitted, among other things to the condition f2 = il exp at p = 
requires to refit the parameters. However the results of Vaks et al 1977 show 
that even in metals as light as Li and Na such refitting is not very essential. 
In addition, in transition metals, in particular in Ir, the choice of the value 
for the effective charge is not quite unumbiguous, and strictly speaking, Z 
can be dependent of (Z = 1 in alkali metals). Because we ignore this 
knowingly more important effect, the results obtained in refitting the pseu- 
dopotential parameters with allowance for p zp are not presented. The direct 
calculation shows that that the change in the results with such a refitting is 
actually small. 

In calculating the major contribution to the pressure po(fl) by (11) in the 
framework of the pseudopotential model for the energy E(Q), use was made of 
the expression derived in the second order of the perturbation theory from a 
pseudopotential with the exchange-correlation energy according to analitical 
approximation (Perdew and Zunger 1981) of the results of the calculation 
by the quantum Montne-Carlo method (Ceperly and Alder 1980). For more 
details see Greenberg et al 1990. The ab initio calculations of the band 
structure and the total energy at different Q were done by the augmentned 
plane wave method. The detailed procedure of these calculations is described 
by Sigalas et al 1992 and Sigalas and Papaconstantopulos 1994. 

§4. Equation of state. Results of calculations 

The strict calculation of equation of state for a metal by Eq. (10) is rather 
difficult. Whereas the electron contributions p and Ap e can be calculated 
form "the first principles" , the calculation of the phonon contributions p zp 
and p p h in the fully "first- principle" approach is a very cumbersome proce- 
dure, and no corresponding results are known to us at that time. For this 
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reason, to calculate the phonon contributions to pressure requires to use some 
model concepts. 

We'll discuss first the results of the T = calculations. Fig. 5 shows the 
values of po(£l) calculated in the framework of the pseudopotential model 
and the band approach that we used. It is seen that the difference of the two 
curves is not too great, at least, at moderate pressures. The reason is that 
the values of the bulk modulus: 3.76 Mbar calculated with the used variant 
of the band theory (Sigalas et al 1992) and 2.93 Mbar in the pseudopotential 
model (see I) do not differ too much. 

Table 2 lists the contribution of the zero point vibrations to the equation 
of state at T = 0. It is seen that this contribution is negligible. The results 
of the calculation of the phonon pressure in the pseudopotential model are 
shown in figs. 6 and 7. The temperature dependent electron contribution to 
the pressure was derived from the results of the band calculations according to 
(13). It appears to be very small in comparison with the phonon contribution 
(50 times smaller at the temperature 295 K). The resulting equation of state 
for Ir at different temperatures is presented in fig. 8. The results for Rh 
are similar and not presented here. Unfortunately, at present we have no 
experimental data to be compared with the above results. 

§5. Elastic moduli at high pressures 

Information about the behavior of the Ir and Rh lattice properties at 
different pressures can be obtained in the framework of the simple pseudopo- 
tential model used here. The extraction of such information from experi- 
ments seems today to be difficult. Therefore, we present the results of the 
corresponding calculations for Ir (fig. 9). The results for Rh are similar. 

§6. Conclusion 

The reasons why the lattice properties of Ir and Rh are so succesfully 
described in the framework of the used simple pseudopotential model are not 
fully understood. However it is believed that if not sufficient, then at least 
necessary condition for this is the absence of any noticable singularities in 
the density of electron states in the immediate vicinity of E F . The results 
of the calculations made by us at different pressures (figs. 10,11) show that 
at not too high pressures this remains the feature of Ir, and, what's more, 
the character of the electron structure near Ep is unchanged. Therefore, it 
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can be assumed that whatever the reason for the succesful pseudopotential 
description at p — 0, the situation at high compressions will not be too 
different from that at p — 0. As a result, one can expect that the information 
about the lattice properies of Ir (and its analogue Rh) at p ^ which was 
derived from the pseudopotential calculations appears ot be raliable. 

The investigation described in the present paper became possible in part 
thanks to the financial support of the International Science Foundation (grant 
RGQ300) and the Russian Foundation for Fundamental Research (grant 95- 
02-06426). 
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List of tables to the paper by Katsnelson et al "Thermal 

EXPANSION ..." 



Table 1. The values of the phonon frequencies cu(q) and the microscopic 
Gruneisen parameters 7(q) in the simmetric points points of the Brillouin 
zone for Ir and Rh. c<j(q) in THz. A-H stands for the values calculated with 
the Animalu-Heine pseudopotential, I stands for the values calculated with 
the addition proposed in I. 

Table 2. The contribution of the zero point vibration to energy and pres- 
sure for Ir and Rh. A-H stands for the values calculated with the Animalu- 
Heine pseudopotential, I stands for the values calculated with the addition 
proposed in I. 
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Figure captions to the paper by Katsnelson et al "Thermal 

EXPANSION ..." 

Fig. 1 Microscopic Gruneisen parameters for Ir calculated for the pseu- 
dopotential from Ivanonv et al 1994 (dashed line) and from I (solid line). 

Fig. 2 Microscopic Gruneisen parameters for Rh calculated for the pseu- 
dopotential from Ivanonv et al 1994 (dashed line) and from I (solid line). 

Fig. 3 Temperature dependence of the linear coefficient of thermal ex- 
pansion (a) and the macroscopic Gruneisen parameter (b) for Ir. Dashed 
line - calculation for the pseudopotential from Ivanov et al 1994; solid line - 
calculation for the pseudopotential from I. o in fig. 3a — experiment from 
KORENOVSKY AND POLYAKOVA, 1990. T pl = huj pi = h(4nZe 2 /Mil ) 1/2 , 
oj p i is the ionic plasma frequency. 

Fig. 4 Temperature dependence of the linear coefficient of thermal ex- 
pansion (a) and the macroscopic Gruneisen parameter (b) for Rh. Dashed 
line — calculation for the pseudopotential from Ivanov et al 1994; solid line 
- calculation for the pseudopotential from I. o in fig. 4a — experiment from 
KORENOVSKY AND Polyakova, 1990. T pl = hu pl = h(4irZe 2 /MQ ) 1/2 , 
ui p i is the ionic plasma frequency. 

Fig. 5 Equation of state for Ir at T = 0. Solid line — calculation for the 
pseudopotential from I; o — ah initio calculations. 

Fig. 6 Dependence of the phonon contribution to pressure upon volume 
at T = 295K and T = T m for Ir (solid line) and Rh (dashed line). The 
calculations were made for the pseudopotential from I. 

Fig. 7 Temperature dependence of the phonon contribution to pressure 
for Ir (solid line) and Rh (dashed line). The calculations were made for 
the pseudopotential from I. T p i = hcu p i = h(4n Z 'e 2 / 'MVLq) 1 ^ 2 , uj p i is the ionic 
plasma frequency. 

Fig. 8 Equation of state for Ir at T=295 K and T = T m . o is ab initio 
calculation of po with pseudopotential calculation of p p h using the addition 
proposed in I at temperature 295 K. o is the same at the temperature T = T m . 
Solid line — all quantities were calculated with the pseudopotential proposed 
in I, temperature is 295 K. Dashed line — the same at the temperature 
T = T 

Fig. 9 Dependence of the elastic moduli Bg~ on volume for Ir. The 
calculations were made for the pseudopotential from I. 
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Fig. 10 The density of electron states (a) and the electron energy spec- 
trum (b) in Ir at equilibrium volume. 

Fig. 11 Density of electron states in Ir at different pressures. 
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table 1 





w(q) 


7(q) 




(1,0,0) 


(1,1/2,0) 


(1/2,1/2,1/2) 


(1,0,0) 


(1,1/2,0) 


(1/2,1/2,1/2) 


Ir A-H 


8.06 


7.00 


8.11 


2.39 


2.27 


2.33 




5.65 


5.85 


4.03 


2.18 


1.83 


1.60 


Ir I 


6.98 


5.90 


7.44 


2.90 


3.01 


2.55 




5.16 


5.22 


4.01 


2.89 


2.17 


1.72 


Rh A-H 


10.03 


8.63 


9.97 


2.17 


2.10 


2.15 




7.04 


7.10 


4.91 


1.99 


1.68 


1.44 


Rh I 


7.62 


6.08 


8.51 


3.12 


3.76 


2.51 




5.93 


5.64 


4.86 


3.34 


2.37 


1.64 



table 2 





Ir A-H 


Ir I 


Rh A-H 


Rh I 


E Z p/E 10 4 


2.49 


2.53 


4.54 


3.78 


p zp /p io 3 


2.63 


3.1 


3.77 


5.25 
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